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We measure the branching fractions and CP asymmetries in the decays B° — » K + K~ K® and 
B + — * K + K° S K% using a sample of approximately 122 million BB pairs collected by the BABAR 
detector. From a time-dependent analysis of the K + K~ K® sample that excludes 4>Ks, the values 
of the CP- violation parameters are S = -0.56 ± 0.25 ± 0.04 and C - -0.10 ± 0.19 ± 0.10, where 
the first uncertainty is statistical, the second is systematic. We confirm that the final state is nearly 
purely CP-even. Using this result and setting C — 0, we extract the Standard Model parameter 
sin 2/3 = 0.57 ± 0.26 ± 0.04t'o' 17 where the last error is due to uncertainty on the CP content. 
We present the first measurement of the CP- violating charge asymmetry Acp(B + — > K + K®K%) = 
— 0.04±0.11±0.02, with a 90% confidence-level interval of [—0.23, 0.15]. The branching fractions are 
B{B° ^K+K~K°) = (23.8±2.0±1.6)xl0~ 6 and B(B+ -^K+K° S K^) = (10.7±1.2±1.0) X 10 -8 . 



PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



In the Standard Model (SM) of particle physics, the de- 
cays B° -> K+K-K° s and B+ -> K+K°K°M are domi- 
nated by b — ► sss gluonic penguin diagrams [2J. CP viola- 
tion in such decays arises from the Cabibbo-Kobayashi- 
Maskawa (CKM) quark- mixing mechanism [3J. Neglect- 
ing CKM-suppressed contributions, the expectation for 
the CP-asymmetry parameters in B° — > K + K~K® de- 
cays is the same as in B° — » JfyK® decays, where CP 
violation has been observed [4J, |5j . The decay rates for 
B+ -> K+K^K° S and B~ -> K~K°K° are expected 
to be equal. However, contributions from physics be- 
yond the SM could invalidate these predictions 6] . Since 
b — * sss decays involve one-loop transitions, they are 
especially sensitive to additional contributions. Present 
results in decays of neutral B mesons are inconclusive 
due to large statistical errors. Belle measures the CP 
asymmetry parameter in <f>Kg decays of sin 2/3 = — 0.96± 
0.50to'?i 7] which is 3.5 standard deviations from the SM 
expectation of sin 2/3 = 0.731 ± 0.056 Jp . A BABAR 
measurement of sin 2/3 = 0.47±0.34+g'Qg |8| is consistent 
with the SM and disagrees with Belle by 2.3 standard 
deviations. 

A more accurate CP measurement can be made us- 
ing all the decays to K + K~K® that do not contain a 
(f> meson. This sample is several times larger than the 
sample of <j>K° HEg. As Belle noted [HJ, the CP con- 
tent of the final state can be extracted using an isospin 
analysis. In decays that exclude <j>K®, Belle measures 
sin 2/3 = 0.51 ± 0.26 ± 0.05±£ 18 0, consistent with the 
SM expectation. In this letter we present measurements 
of CP asymmetry and CP content in K + K~ K® decays, 
and the first measurement of the charge asymmetry rate 
in B+ -y K+K° s K a s decays. 

This analysis is based on about 122 million BB pairs 
collected with the BABAR detector [ll| at the PEP-II 
asymmetric-energy e + e~ storage rings at SLAC, operat- 
ing on the T(AS) resonance. We reconstruct B mesons 
from K® — ► 7r + 7r - and K ± candidates. Charged kaons 
are distinguished from pions and protons using energy- 
loss (dE/dx) information in the tracking system and from 
the Cherenkov angle and number of photons measured 
by the detector of internally reflected Cherenkov light 
(DIRC). We accept K° — ► ir + ir~ candidates that have a 
two-pion invariant mass within 12 MeV/c 2 of the nomi- 



nal K^ mass [T2J , a decay length greater than 3 standard 
deviations, and a cosine of the angle between the line 
connecting the B and K® decay vertices and the K® mo- 
mentum greater than 0.999. The three daughters in the 
B decay are fitted constraining their paths to a common 
vertex, and the K® mass to the nominal value. 

In the characterization of the B candidates we use 
two kinematic variables. The energy difference AE = 
Eb — a/s/2 is reconstructed from the energy of the B can- 
didate Eb and the total energy yfs in the e + e _ center-of- 
mass (CM) frame. The AE resolution for signal events is 
18 MeV. We also use the beam-energy-substituted mass 
toes = \/( s / 2 +^ -PB) 2 /Ef -p|, where ($,£*) is the 
four-momentum of the initial e + e~ system and pb is the 
momentum of the B candidate, both measured in the 
laboratory frame. The toes resolution for signal events is 
2.6 MeV/c 2 . We retain candidates with \AE\ < 200 MeV 
and 5.2 < to es < 5.3 GeV/c 2 . 

The background is dominated by random combinations 
of tracks created in e + e~ — ► qq (q = u, d, s, c) continuum 
events. We suppress this background by utilizing the dif- 
ference in the topology in the CM frame between jet-like 
qq events and spherical signal events. The topology is de- 
scribed using angle 9t between the thrust axis of the B 
candidate and the thrust axis of the charged and neutral 
particles in the rest of the event (ROE) [111 . Other quan- 
tities that characterize the event topology are two sums 
over the ROE: L = J2\p">*\ an d L 2 = J2 \Pi*\ cos 2 9i, 
where #.; is the angle between the momentum pi* and 
the thrust axis of the B candidate. Additional sepa- 
ration is achieved using the angle 9b between the B- 
momentum direction and the beam axis. After requiring 
| cos6>t| < 0.9, these four event shape variables are com- 
bined into a Fisher discriminant T [13J. 

The remaining background originates from B decays 
where a neutral or charged pion is missed during recon- 
struction (peaking B background). We use Monte Carlo 
(MC) events to model the signal and the peaking back- 
ground, and data sidebands to model continuum back- 
ground. 

We suppress background from B decays that proceed 
through a b — ► c transition leading to the K + K~K® 
(K + K®K®) final state by applying invariant mass cuts 
to remove D° -> KK, D+ -> K+K°, J/ip -* KK, and 



XcQ — * KK(K^K^) decays. Finally, B decays into final 
states with pions are eliminated by requiring the pion 
misidentification rate to be less than 2%. 

The time-dependent CP asymmetry is obtained by 
measuring the proper time difference At between a fully 
reconstructed neutral B meson (Bcp) decaying into 
K + K~ K®, and the partially reconstructed recoil B me- 
son (Stag)- Decay products of the recoil side are used 
to determine the B tag meson's flavor (flavor tag) and to 
classify the event into five mutually exclusive tagging cat- 
egories 4]. If the fraction of events in category c is e c 
and the mistag probability is w c , the overall quality of 
the tagging, £ c e c (l - 2w c ) 2 , is (28.0 ±0.4)%. 

The time difference At is extracted from the measure- 
ment of the separation Az between the Bqp and B tag 
vertices, along the boost axis (z) of the BB system. The 
vertex position of the Bqp meson is reconstructed pri- 
marily from kaon tracks, and its MC-estimated resolu- 
tion ranges between 40-80 /im, depending on the open- 
ing angle and direction of the kaon pair. The final At 
resolution is dominated by the uncertainty on the -Btag 
vertex which allows the At (Az) precision with r.m.s. of 
1.1 ps (180 /im). We retain events that have |Ai| < 20 ps 
and whose estimated uncertainty ua* is less than 2.5 ps. 
The At resolution function is parameterized as a sum 
of two Gaussian distributions whose widths are given by 
a scale factor times the event-by-event uncertainty a At- 
A third Gaussian distribution, with a fixed large width, 
accounts for a small fraction of outlying events |4| . 

Parameters describing the tagging performance and 
the At resolution function are extracted from ap- 
proximately 30,000 B° decays into D^~X+ (X+ = 



■ P" 



flavor eigenstates (-Bfl av sample). 



The decay rate f + (f _ ) when the flavor of the tagging 
meson is a B° (B°) is given by 



f±(At) 



-|At|/r B o 



4t b o 

T C cos ( Am d At)} 



{ l±Ssm(Am d At) 



(1) 



where t b o is the mean B° lifetime and Amj is the B°- 
B° oscillation frequency. The parameters C and S de- 
scribe the magnitude of CP violation in the decay and 
the interference between decay and mixing, respectively. 
In the SM, we expect C — because there can be no 
direct CP violation when there is only one decay mech- 
anism. If we exclude 4>K® events by applying a K + K~ 
invariant mass cut of 15 MeV/c 2 around the nominal <fi 
mass |l2| . and assume that the remaining Bcp candi- 
dates are CP-even, as our analysis below indicates, we 
expect S= -sin 2/3 = -0.731 ±0.056 0,0. 

Direct CP violation in B + — > K + K®K® decays is mea- 
sured as an asymmetry in the decay rates 



*-K~K%K% 

•Anp = 



K-K%K% 



K+K"K" 



(2) 



The SM expectation for Acp is zero. 

Branching fractions and CP asymmetries are extracted 
in unbinned extended maximum likelihood fits to the dif- 
ferent samples. The likelihood function C, with event 
yields iVj and probability density functions (PDFs) Vij, 
is: 



£=«p -5>i n 



£ww 



(3) 



where j runs over events and i over event yields. Wc 
have a total of 6144 events in the K+K2.K® mode, and 
13864 (12862) in the K + K~K° mode with (j)K° in- 
cluded (excluded). 

In the measurement of the branching fractions B, the 
total PDF is formed as V{m ES ) • V{AE) ■ V(T). Event 
yields for signal, continuum, and peaking B background 
are varied in the fit. In the extraction of the charge asym- 
metry Acp in K + K®K® decays, the yields are split by 
the charge, which brings the total number of varied pa- 
rameters to six. To extract the branching fractions, we 
assign a weight for each event to belong to the signal 
decay, Wj — y ^-'.p.'' J where V S) i is the signal row of 
the covariance matrix obtained from the fit 14]. The 
branching fraction is calculated as B = J2j Wj/ej. Since 
the efficiency £j varies across the phase space, Ej is com- 
puted in small phase-space bins using simulated events. 
The method is cross-checked with a simple counting anal- 
ysis. Distributions of toes and AE are shown in Fig. ^ 
and the fit results are given in Table [I] 
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FIG. 1: Projection plots of the variables mES (a, c) and AE 
(b, d) in the fits for B° -> K + ' K~ K° s (top) and B + -> 
K + K®Ks (bottom) decays. The points are data and the 
curves are projections from the likelihood fit. The signal-to- 
background ratio is enhanced with a cut on the event proba- 
bility. 

In the time-dependent CP fit, K + K~K® events that 
exclude (j>K® decays are fit simultaneously with the -Bfl av 
sample. The PDFs are formed as V(m ES )-V(AE)-'P(J r )- 
V c (At; a At) for B C p events and V(m ES ) ■ Pc(At; a At) in 



TABLE I: Summary of branching fraction (B), time-dependent (S, C) and direct CP-asymmetry (Acp) results. iV s ig and e are 
the signal yield and the average total efficiency in the branching-fraction fit; f eV en is the CP-even fraction of the final states. 
The 90% confidence-level interval for Acp is [-0.23,0.15]. 



Mode 


e(%) 


^ig 


B (lO" 6 ) 


J even 


S 


C 


A C p 


K+K-K CP 
K+K-K a11 


8.58 
8.78 


201 ± 16 
249 ± 20 


20.2 ±1.9 ±1.4 
23.8 ±2.0 ±1.6 


0.98 ±0.15 ±0.04 
0.83 ±0.12 ±0.03 


-0.56 ±0.25 ±0.04 


-0.10 ±0.19 ±0.10 


— 


K+K°K° 


9.7 


122 ± 14 


10.7 ±1.2 ±1.0 


— 


-0.16 ±0.35 


-0.08 ±0.22 


-0.04 ±0.11 ±0.02 



CP Excludes <j>K° s events. 



the -Bfiav sample. The At resolution and tagging param- 
eters are allowed to be different for each tagging category 
c. Fit parameters that are common to both samples are 
the signal fractions in tagging categories e c , the average 
mistag fraction w c , the difference between B° and B° 
mistag rates Aw c , and the At resolution functions for sig- 
nal and background events. We also vary the K + K~ K® 
signal yield and background yields in tag categories, the 
CP parameters, and the parameters describing the At 
shape of the background. The total number of floated 
parameters is 38. The largest correlation between S or 
C with any linear combination of other parameters is 
6.6%. 




-6-4-20246 
At ( P S) 

FIG. 2: Plots a) and b) show the At distributions of B°- and 
_B°-tagged K + K~ K® events. The solid lines refer to the fit 
for all events; the dashed lines correspond to the background. 
Plot c) shows the raw asymmetry, where the solid line is ob- 
tained from the fit and the dotted line is the SM expectation 
for the measured CP content. The signal-to-background ratio 
is enhanced with a cut on the event probability. 

Results of the time-dependent CP asymmetry mea- 
surement in K + K~K® are given in Table |U Figure El 
shows the At distributions of events with B° and B° 
tags, with projections from the likelihood fit superim- 
posed. The fit procedure is verified with the K + K®K® 
sample (TableHJ, where one expects zero asymmetry, and 
the J/ipK® sample where the results are consistent with 
our previous measurement [4|. 

We evaluate the fraction / eV on of CP-even final states 
in B° — > K + K~K® decays by comparing K + K~K° and 

K+K s K s decay rates: /even = 2 T(B«^K+K-K») E3' 



The results listed in Tableware in agreement with Belle's 
measurements of 0.86 ±0.15 ±0.05 and 1.04 ±0.19 ±0.06 
for the total sample and the CP sample that excludes 
(j)K® events, respectively |l(J- We estimate the fraction 
of remaining 4>K® events in the CP sample, using a non- 
interfering Breit-Wigner for the </> shape and measured 
branching fractions, to be 1.1 ± 0.4%. As a consistency 
check, we examine the distribution of the cosine of the 
helicity angle 9h, which is defined as the angle between 
the K + and B° directions in the K + K~ center of mass 
frame. The distribution in several K + K~ invariant mass 
bins of the CP sample is approximately uniform which 
is consistent with S-wave decays. The presence of in- 
terference effects due to CP-odd amplitudes cannot be 
ruled out, but this would require a full amplitude analy- 
sis which is not feasible with the present statistics. 

If we account for a small CP-odd fraction in the CP 
sample, we can extract the SM parameter sin 2/3. In a fit 
with C=0we get sin2/3 = -S/(2f even - 1) = 0.57 ± 
0.26±0.04^q' 17 where the last error is due to uncertainty 
on the CP content. 



TABLE II: Branching fraction systematic uncertainties (%). 



Source 



K+K-K". K + K%K 



TW^O" 



Efficiency 


5.6 


8.6 


PDF parameterization 


2.7 


2.5 


Non-charm BB background 


2.2 


2.9 


Charm BB background 


1.2 


1.0 


Other 


1.7 


1.6 



Total 



6.9 



0.6 



Systematic uncertainties in the branching fraction 
measurements are given in Tabled We include contribu- 
tions from the signal reconstruction efficiency and from 
the modeling of the efficiency variation over the phase 
space. Other errors come from the fit bias, the count- 
ing of BB pairs, and the misidentification of kaons. We 
assume equal production rates of B°B° and B + B~ . The 
systematic uncertainty on Acp due to charge asymmetry 
in track finding and identification is 0.02. 

The systematic errors on the time-dependent CP- 
asymmetry parameters are given in Table lllll The errors 
account for the fit bias, the presence of double CKM- 
suppressed decays(DCSD) in P tag [l3, uncertainty in the 
beam spot and detector alignment, and the asymmetry in 



TABLE III: Systematic uncertainties in CP parameters. 



boldt Foundation. 



Source 



C 



Fit bias 0.024 0.026 

DCSD 0.018 0.053 

Detector effects 0.013 0.012 
Tag asymmetries 0.010 0.078 
Other 0.016 0.012 



Total 



0.04 0.10 



the tagging efficiency for signal and background events. 
Other smaller effects come from At resolution, PDF pa- 
rameterization of yield variables, and uncertainty on the 
B° lifetime and mixing frequency In the fit we use Tgo = 
1.537 ± 0.015 ps and Am d = 0.502 ± 0.007 ps' 1 0. 

In summary, we have measured branching fractions for 
charmless decays of B mesons into the three-body final 
states B° -> R+R-R and B+ -► R+R°R° S . Using 
two independent approaches, we find that the R + R~ R® 
final state is dominated by a CP-even component. The 
results agree with previous measurements yj |lj| . In the 
first measurement of the charge asymmetry in B + — v 
R + R°R° decays, we find no evidence for direct CP vi- 
olation. We measure a time-dependent CP asymmetry 
in B° -> R+R-R° s decays at the 1.9a level. The ob- 
tained sin 2/3 is consistent with the SM expectation and 
previous measurements in decays into the R + R~R® fi- 
nal state [Ilia, but differs from Belle's measurement in 
4>R° decays Q by 2.7 standard deviations. 

We are grateful for the excellent luminosity and ma- 
chine conditions provided by our PEP-II colleagues, and 
for the substantial dedicated effort from the computing 
organizations that support BABAR. The collaborating 
institutions wish to thank SLAC for its support and 
kind hospitality. This work is supported by DOE and 
NSF (USA), NSERC (Canada), IHEP (China), CEA and 
CNRS-IN2P3 (France), BMBF and DFG (Germany), 
INFN (Italy), FOM (The Netherlands), NFR (Norway), 
MIST (Russia), and PPARC (United Kingdom). Indi- 
viduals have received support from the A. P. Sloan Foun- 
dation, Research Corporation, and Alexander von Hum- 



Now at Department of Physics, University of Warwick, 
Coventry, United Kingdom 
t Also with Universita della Basilicata, Potenza, Italy 
* Also with IFIC, Instituto de Ffsica Corpuscular, CSIC- 

Universidad de Valencia, Valencia, Spain 
s Deceased 
[1] Charge-conjugate states are included unless explicitly 

stated otherwise. 
[2] D. London and R. D. Peccei, Phys. Lett. B 223, 257 
(1989); N. G. Deshpande and J. Trampetic, Phys. Rev. 
D 41, 895 (1990); R. Fleischer, Z. Phys. C 62, 81 
(1994); N. G. Deshpande and X. G He, Phys. Lett. B 
336, 471 (1994); Y. Grossman, Z. Ligeti, Y. Nir and 
H. Quinn, Phys. Rev. D 68 (2003) 015004 ; M. Gronau 
and J. L. Rosner, Phys. Lett. B 564, 90 (2003). 
[3] N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963); M. Ko- 
bayashi and T. Maskawa, Prog. Theor. Phys. 49, 652 
(1973). 
[4] B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 

89, 201802 (2002). ^_^ 

[5] K. Abe et al. [Belle Collaboration], hep-ex/0308036 

[6] Y. Grossman and M. P. Worah, Phys. Lett. B 395, 241 

(1997); R. Fleischer, Int. J. Mod. Phys. A 12, 2459 

(1997); D. London and A. Soni, Phys. Lett. B 407, 61 

(1997). 

[7] K. Abe et al. [Belle Collaboration], Phys. Rev. Lett. 91, 

261602 (2003). ' _^ 

[8] B. Aubert [BABAR Collaboration], hep-ex/0403026 
[9] B. Aubert et al. [BABAR Collaboration], Phys. Rev. D 
69, 011102 (2004). 
[10] A. Garmash et al. [Belle Collaboration], Phys. Rev. D 

69, 012001 (2004). 
[11] B. Aubert et al. [BABAR Collaboration], Nucl. Instrum. 

Meth. A 479, 1 (2002). 
[12] K. Hagiwara et al. [Particle Data Group Collaboration], 

Phys. Rev. D 66, 010001 (2002). 
[13] R. A. Fisher, Annals Eugen. 7 179 (1936). 
[14] M. Pivk, F. R. Le Diberder, physics/0402083 
[15] O. Long, M. Baak, R. N. Cahn and D. Kirkby, Phys. 
Rev. D 68, 034010 (2003). 



We measure the branching fractions and CP asymmetries in the decays B° — > K + K~ K$ 
and B + — > K + K$Ks using a sample of approximately 122 million BB pairs collected by 
the BaBar detector. From a time-dependent analysis of the K + K~K$ sample that excludes 
4>Kg, the values of the CP-violation parameters are S = —0.56 ± 0.25 ± 0.04 and C = 
—0.10 ± 0.19 ± 0.10, where the first uncertainty is statistical, the second is systematic. We 
confirm that the final state is nearly purely CP -even. Using this result and setting C = 0, 
we extract the Standard Model parameter sin 2/5 = 0.57 ± 0.26 ± 0.04^q' 17 where the last 
error is due to uncertainty on the CP content. We present the first measurement of the 
CP- violating charge asymmetry A C v(B + -► K + K S K S ) = -0.04 ± 0.11 ± 0.02, with a 90% 
confidence-level interval of [—0.23,0.15]. The branching fractions are B(B° — > K + K~ K°) = 
(23.8 ± 2.0 ± 1.6) x 10~ 6 and B(B + -* K + K S K S ) = (10.7 ± 1.2 ± 1.0) x 10~ 6 . 
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